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ABSTRACT 

We present results obtained from a set of cosmological hydrodynamic simulations of galaxy clusters, 
aimed at comparing predictions with observational data on the diversity between cool-core (CC) and non¬ 
cool-core (NCC) clusters. Our simulations include the effects of stellar and AGN feedback and are based 
on an improved version of the smoothed particle hydrodynamics code GADGET-3, which ameliorates gas 
mixing and better captures gas-dynamical instabilities by including a suitable artificial thermal diffusion. 
In this Letter, we focus our analysis on the entropy profiles, the primary diagnostic we used to classify the 
degree of cool-coreness of clusters, and on the iron profiles. In keeping with observations, our simulated 
clusters display a variety of behaviors in entropy profiles: they range from steadily decreasing profiles 
at small radii, characteristic of cool-core systems, to nearly flat core isentropic profiles, characteristic of 
non-cool-core systems. Using observational criteria to distinguish between the two classes of objects, we 
find that they occur in similar proportions in both simulations and in observations. Furthermore, we 
also find that simulated cool-core clusters have profiles of iron abundance that are steeper than those 
of NGG clusters, which is also in agreement with observational results. We show that the capability 
of our simulations to generate a realistic cool-core structure in the cluster population is due to AGN 
feedback and artificial thermal diffusion: their combined action allows us to naturally distribute the 
energy extracted from super-massive black holes and to compensate for the radiative losses of low-entropy 
gas with short cooling time residing in the cluster core. 

Subject headings: galaxies: clusters: general - galaxies: clusters: intracluster medium - X-rays: 
galaxies: clusters - methods: numerical 


1. INTRODUCTION 

During their hierarchical assembly, galaxy clusters grow 
in mass via diffuse accretion processes as well as via merger 
events. These processes lead to the shock heating of the 
intra-cluster medium (IGM) to a virial temperature of up 
to T ~ 10® K, with central electro n number d ensities cor- 

Wit|2005| ^ 
is density, hi 


responding to n, 


10 ®cm ® 


Borgani 2012 


(e-g- 


Kravtsov & 


rot baryons 


__ for a review). At this density, 

in the core regions have a cooling time, Gooi cx; T^/^/rie, 
shorter than the Hubble time. Under this condition, the 
ICM suffers from radiative losses and becomes colder and 
denser. This process might initiate a run-away cooling 
process, which would cause an extreme accretion rate in 
the central galaxy, unless some heating mechanism were 
to balance the cooling. Substantial evidences from X-ray 
observations proved that the regulating mech anism is the 
feedback by active galactic nuclei (AGNs, e.g. |McNamara 
fc Nulsen||2007 Voit et aL]|2015 ). 


TE a situation of steady accretion, low-entropj0and metal- 
enriched gas associated with merging substructures, or 
tunneled by filaments, sink toward the central regions of 
the so-called cool-core clusters (GG, named as such by 


^ We adopt the standard in the X—ray definition of entropy: K = 
Twith T being the ICM temperature and rie being the electron 
number density. 


Molendi & Pizzolato 2001|). On the contrary, non-cool- 
core (NCC) objects are observed with nearly isentropic 
gas cores at a higher entropy level and show no evidence 
of a spike in their rnetal a bundances (|De Grandi et al. 
Maughan et al.||2008 ). 


2004 


While idealized simulations reproduce some of the many 


featu res of the CG systems (e.g Gaspari 2015 Li et al. 


20151, producing the diverse populations of CC and NCC 


dusters has so far proved to be so far quite a formidable 


challenge for cosmologic al hydrodynamic simulations (Bor¬ 
gani & Kravtsov 2011 and references therein). Earlier 
simulation works predicted that the origin of the two classes 
of CC and NCC clusters is “primeval”, with CC being 
destroye d exclusively by early mergers or preheating phe¬ 
nomena ( B urns et al. 2008 McCarthy et al. 2008 jPoole 


et al. 2008 Planelles & Quilis 2009|). In fact, a CC cius- 
ter was claimed to maintain its high-density and metal- 
rich core even after a major merger, if that happens at 
z < 0.5. However, this prediction is in conflict with ob¬ 
servations and recent cosmological simulations by|Hahn~et 


al. (20151. In most cases, CC clusters show regular X-ray 
morphology while the opposite is true for NCC systems. 
Data, thus, favor an “evolutionary” model in which a tran¬ 
sition between the CC and NCC status of a cluster can 
take place over relatively short timescales in consequence 
of a merger event that mixes the convective stable gas and 


1 
































































2 


Rasia et al. 


destroys the CC (Rossetti et al. 2011). Specifically, Hahn 
et al. (2015) advocate that the dichotomy between CC and 
NCC clusters originates by low-angular-momentum major 
mergers. 

In this Letter, we present a set of simulated clusters in 
which, for the first time, we find (1) co-existence of CC 
and NCC systems, (2) agreement with the observed ther¬ 
modynamical and chemical properties of the two classes of 
clusters, and (3) evidence of a recent transition between 
CC and NCC in the evolution of individual objects. As we 
will discuss, this result is achieved thanks to the combined 


action of a new AGN feedback model (Steinborn et al. 
2015) and an i mproved smoothe d-particle-hydrodynamics 
(SPH) scheme ( Beck et ar]|2015 ). 


2. SIMULATIONS 

We provide a short description of the simulations ana¬ 
lyzed here, while we refer to a forthcoming paper (Planelles 
et al. 2015, in preparation) for a more thorough descrip¬ 
tion. The simulations have been generated from the same 
set of initia l conditions originally presented by [Bonafed^ 
et al. (2011). They correspond to zoomed-in initial condi¬ 
tions of 29 Lagrangian regions selected around 24 massive 
clusters with M 500 between 5 and 20 xlO^"^ h~^MQ and 5 
poorer clusters with M 500 in the range 0.7-3xl0^‘^h~^M q. 
The simulations are carried out with the GADGET-3 code 


(Springel 2005) and include an upgraded versio n of the 
SPH scheme, as described in I Beck et al. (2015). These 
developments consist of (a) a higher-order (Wendland C'^ 
instead of the standard B-spline) interpolating kernel to 
better describe discontinuities and reduce dumpiness in¬ 
stability, (b) a time-dependent artificial viscosity term to 
capture shocks and minimize viscosity away from shock 
regions, (c) an artificial conduction term that largely im¬ 
proves the SPH capability of following gas-dynamical in¬ 
stabilities and mixing processes. With respect to the phys¬ 
ical thermal conduction, it similarly promotes the trans¬ 
port of heat, since it is applied to contact discontinuities 
in internal energy, but, differently, it does not depend on 
thermodynamical properties of the gas. To assess the per¬ 
formances of this improved SPH implementation, several 
standard hydrodynamical tests, such as weak and strong 
shocks, shear flows, gas mixing, and self- gravitating gas, 
have been performed in Beck et al. (2015). 

The physical processes included in our simulations can 
be described as follows. Metallicity-dependent radiative 
cooling and the effect of a uniform time-dependent UV 
backgro und are considered a s in Planelles et al. (2015, 
see also Wiersma et al.|[2009 ). A sub-resolution model for 
star formation fro m a multi-phase interstellar m edium is 
implemented as in Springel & Hernquist (2003). Kinetic 
feedback driven by supernovae (SN) is accounted in the 
form of galactic winds with a velocity of ~ 350 km s“^. 
Metal production from SN-H, SN-Ia and asymptotic-giant- 
branch (AGB) stars follo ws the original receip t by Torna- 
tore et al. (2007, see also |Planelles et al.|2014 ). The AGN 
feedback is modeled with the implementations recently 


presented in Steinborn et al. (2015). The released ther¬ 
mal energy accounts for contributions by both mechanical 
outflows and radiation, separately computed in the code. 
Eddington-limited gas accretion onto BHs is computed by 
multiplying the Bondi rate by a boost factor a. In this 


Letter, we consider only cold accretion (ttcoid = 100 and 
cihot = 0). We verified that the results hold when hot 
accretion is included with ahot = 10. For the feedback en¬ 
ergy, we assume variable efficiencies for the radiation and 


the mechan ical outflow (Equations 19 and 20 of Steinborn 
et al. 2015) that depend on the BH mass and on the ac- 


cretion rate. Finally, we fix e/ = 0.05 as the efficiency of 
coupling the radiated BH energy to the IGM. 

We emphasize that the parameters defining the modi¬ 
fication of SPH and the AGN feedback model have been 
tuned to overcome the limitations of standard SPH in pass¬ 
ing classic hydrodynamical tests and to reproduce the ob¬ 
servational scaling relation between BH mass and stellar 
mass of t he host galaxies, includ ing the brightest cluster 
galaxies (McConnell & Ma 2013). No attempt has been 
pursued to reproduce any of the observational properties 
of the ICM. The existing match between other simulated 
and observed stellar and hot gas properties will be pre¬ 
sented in future papers. 

The dark matter mass resolution is equal to 8.3x 10 ®/i“^Mq 
and the z = 0 Plummer-equivalent softening length for the 
gravitational force is e = 3.75h~^ kpc. In our sample, this 
value is always inferior than 0.01 x R^qq. 

3. RESULTS 

Several approaches have been used in the literature to 
classify clusters as CC and NCC systems based, e.g., on the 
central gradient of temperature or entropy, on the central 
entropy level, on the cooling time, on the mass deposition 
rate inferred from X-ray spectroscopy, on the cuspiness of 
the central gas density profile (comparisons among these 
criteria are presented in [Cavagnolo et al.||200^ |Hudson| 
et al.||2010[ [Leccardi et ^ |2010[ [McDonald et al.||2013| 
Pascut & Ponman 2015). In order to assign a degree of 
“cool-coreness” to our simulated clusters, we jointly use 
the following two criteria which measure the shape and 
level of the entropy profiles in the central regions and are 
extensively applied to observation al data. 

Following Rossetti et al. (2011, see also Leccardi et al. 
2010 ), we base the first criterion on the computation of 
the pseudo entropy: 

^ (^in/^out) 

(EMin/EMout)^^^ ’ 


where the spectroscopic-like temperature (Mazzotta et al. 
2004), T, and the emission measures, EM, are computed 


in the IN region, r < 0.05 x i?i 8 o, and in the OUT region, 
0.05 X i?i 8 o < r < 0.15 X i?i 8 o- We define CC those clusters 
with a < 0.55. 

In addition, we request the value of the central entropy 
to be Kq < 60 keV cm^. We model the entropy profile as 

K = Ko + K,oo>^{R/Riookr, ( 2 ) 

where we treat as free parameters the core entropy, Kq, 
the entropy at Riook(= 100 kpc), A^ioo, and the slope, 
a, (Cavagnolo et al. 2009). For each simulated cluster the 
entropy prohle is measured in linearly equispaced spherical 
shells, with the innermost radius being that of the sphere 
containing at least 100 gas particles. In all cases, this 
radius is greater than 5 x e. We use Rqqq as the outermost 
radius out to which we perform the fit. Three clusters 
with a < 0.55 have steadily declining entropy profiles. In 
these cases, Kq is not constrained. We, thus, consider the 
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Fig. 1.— Si mulated 2 = 0 entro py (top panels) and Iron abundance (bottom panels) profiles compared with observations by|Pratt et al.| 
l |2010[ | and by |Ettori et al.| | |2015[ l, respectively. For the simulated CC (left panels) and NCC (right panels) profiles: the medians are m 
black and the and the 84*'* percentiles are deli mited by shaded regions. The observed profiles are in cyan (CC) and magenta (NCC). 
Metallicity profiles are expressed in solar units using [Anders &: Greves^||1989|. 


overall profiles and keep in the CC sample the two objects 
with a ~ 1 and exclude that with a < 0.5. 

The distribution of cr and Kq is continuous and no clear 
bimodality i s pres ent, in contrast with recent results by 
Hahn et al. (2015). Considering the aforementioned ob¬ 
servational thresnolds, we find that 38% (11/29) of our 
simulated clusters at z = 0 are classified as CC. This 
fraction is quite close to the 35% value reported by |Eck- 


ert et al. (2011) from a similar mass distributed sample, 
HlCl'LUGS, after accounting for the sample incomplete¬ 
ness and the X-ray selection bias. 

In the upper panels of Figure we compare the en¬ 
tropy profiles of ou r simulated cluster s with observational 
results obtained by Pratt et al. (|2010|) from XMM-Newton 
observations of the KFXCFSS clusters. Our simulations 
now produce a population of CC clusters, whose entropy 
profiles decline down to the innermost regions. This result 
confirms that the energy extracted from the central BHs is 
distributed in such a way to compensate for the radiative 
losses of gas at low entropy and, therefore, to keep it in 
the hot phase despite the fact that it formally had a short 
cooling time. More generally, we note that observed and 
simulated entropy profiles agree quite well in normaliza- 
tion and slope, for both CC and NCC populations. 


of a Virgo-like cluster by including the effect of AGN feed¬ 
back. However, in the same paper, these authors claimed 
that this conclusion is spoiled as soon as the effect of metal¬ 
licity is included in the cooling function. This confirms the 
fragility of a self-regulated AGN feedback. We, therefore, 
compare the simulated mass-weighted iron profiles to the 


metallicity of nearby clusters, z < 0.2, from 

Ettori et al. 

( 

2015P in the bottom panels of Figure|l| In 

seeping with 


ized by a rather high enrichment level in the central regions 
while NCC objects have a shallower peak. This agreement 
suggests that our simulations provide an effective descrip¬ 
tion of the mixing related to gas-dynamical processes act¬ 
ing during the cluster formation and of the rise of enriched 
gas caused by feedback processes: the same processes lead¬ 
ing to the diversity of entropy profiles are also responsible 
for the diversity of the profiles of metal abundance. 

From an object-by-object investigation of our simulated 
clusters, we recognize a variety of situations for the estab¬ 
lishment of a CC and for the transformation between CC 
and NCC, which can take place over a quite a short time 


Dubois et al. (2011) were also able to produce the cor- 
rect CC entropy structure in a high-resolution simulation 


^ In their classification of CC and NCC, Ettori et al. 12015 I used a 
different definition of the IN and OUT regions m the pseudoAntropy 
calculation. For their sample, we consider the limit cr < 0. 6 that 
corresponds to our adopted threshold (cr < 0.55) according tojBaldij 

jet al.jpoT^. 
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31 32 34 38 48 65 101 173 316 


Fig. 2.— Maps of pseudo entropy of two simulated clusters with masses = 2.4 X (upper panels) and M500 = 7.3 X 

10'^^/r“^MQ(lower panels) at z = 0. The size of the images is 1 Mpc and the line-of-sight integration is for 10 Mpc. 


and even at low redshift. As an example, in Figure we 
show the evolutionary sequences of pseudo-entropy maps 
from 2: = 0.8 to z = 0 for two objects. Following an ob¬ 
servational approach, we create maps of spectroscopic-like 
temperature, TWtslj and of X-ray soft band ([0.5-2] keV) 

emissivity, A4x, which we then combine as A^tsl/A^V ^ 
to build maps of pseudo entropy ( Finoguenov et al.|2010| . 
The upper panel shows the case of an NCC cluster at 
z = 0.8 that slowly and continuously develops a CC struc¬ 
ture. On the contrary, the lower panel corresponds to a 
system with a low-entropy core at high redshift. The core 
is still in place at z = 0.25 and it is eventually destroyed 
by a merger between the main halo and a sub-clump at 
z ^ 0.1. Remarkably, this last case of CC remnant is not 
an isolated event: almost half of NCC clusters at z = 0 
were CC at z = 1 ( Rossetti et al.|[20TT ). 

These results vary somewhat from previous claims, also 
based on simulations, which stated that CCs can only be 
destroyed at high redshift through major mergers, while 
being much more resilient against mergers at low redshift 


Burns et al. 2008 


Poole et al. 2008). A close com- 


le.g. _ . _ 

parison of our results with these previous analyses is not 
straightforward due to substantial differences in both the 
hydro dynamic schemes and the implementation of feed¬ 
back. Poole et al. (2008) carried out non-cosmological sim- 
ulations of isolated two-body merger s using a mor e stan- 
dard implementation of SPH, while Burns et al. (2008) 


used a Eulerian code with a cosmological set-up. Neither 
one of these two analyses included the effect of AGN feed¬ 
back and, for this reason, the strength of their cores is 
probably affected by overcoo ling. In order to surpass this 
limitation. Burns et al. (2008 ) truncated star form ation be¬ 
low a certain redshift, while Poole et al. (2008) excluded 
the central 40 kpc in their analysis. Both strategies, how¬ 
ever, do not directly attack the main problem of the del¬ 
icate balance between cooling and heating. This balance, 
essential for the creation of a low-entropy core, needs to 
be regulated since early times and can be achieved in our 
simulations only by including AGN feedback. 

4. DISCUSSION AND CONCLUSIONS 

It has long been known that ICM evolution purely driven 
by gravitational processes leads to self-similar entropy pro¬ 
files that scale as a power-law of radius, K{r) oc r°‘ with 
1 -1.1 i n the cluster outskirts ( Voit et ar]|2005 Nagai 
et al. [2007 ). Non-radiative simulations, however, differ in 
their behavior in the centra l regions. A flat co re profile is 
found in grid-based codes ( Frenk et al. 1999) and in re¬ 
cent improved versions of SPH (Biffi & Valdarnini 2015 


Sembolini et al. 2015 and references therein). Standarc. 


SPH simulations, instead, produce entropy profiles that 
are declining toward the center. 

In addition to mixing, radiative processes have the coun¬ 
terintuitive effect of also raising entropy in core regions 
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(e.g. Borgani & Kravtsov 2011). If not counteracted by 
some energy feedback process, radiative cooling selectively 
removes gas with low entropy and short cooling time from 
the hot-X-ray emitting phase, thus fuelling star formation 
(e.g. Li & Bryan 2012) and leaving in the ICM only gas 
with a relatively higher entropy. 

The AGN feedback plays a key role in explaining our 
results. Considering the new SPH set-up but excluding 
AGN feedback, we find that overcooling leads to a thermal 
structure of central cluster regions whose entropy profiles 
behave as in NCC clusters (top panel of Figure 3). At the 
same time, the excess of star formation produces an ex¬ 
ceedingly high level of ICM metal enrichment, with heavy 
elements mostly locked around galaxies. As a consequence, 
simulated clusters without AGNs always display spikes of 
central metallicity, an enrichment pattern typical of CC 
clusters (bottom panel of Figure 3). As observations sug¬ 
gest, the AGN feedback should have the threefold effect of 
(1) compensating radiative losses so as to prevent gas with 
a short cooling time from dropping out of the ICM, (2) re¬ 
ducing star formation in massive halos at low reds hift, (3) 
raising metals away from star-forming regions (e.g. ^jacki 


et al.|[2007 Fabjan et al.||201~ McCarthy et al.|[20iip. 

However, attempts to create realistic thermal structures 

of CC in simulated clusters so far have failed, even when 
including AGN feedback, which is otherwise successful in 


reproducing X-ray scaling relations (e.g. Puchwein et al. 


2008 
et ai. 


Planelles et al. 2014) and their scatter (Le Brun 


2014p. 'The encouraging success of the siinuiations 


presented here in producing thermal and chemo-dynamical 
structure of CC in agreement with observations is related 
to the combined action of the AGN feedback model and of 
the artificial thermal diffusion. 

Figure |4 compares the entropy prohles of the cluster 
shown in the upper panels of Figure with those of the 
same object simulated without artificial conduction but in¬ 
cludes other recent implementations such as the new kernel 
and the artihcial viscosity. By excluding artificial conduc¬ 
tion, the entropy profiles at z = 0.8 and z = 0 are similar 
to each other with a large isentropic core, characteristics 
of a strongly NCC object. This contrast with the previous 
hnding showing a transition from an NCC to a CC situa¬ 
tion (Figure^. Besides confirming the results reported in 


previous works based on standard SPH (Sijacki et al. 2008 
Planelles et al. 2014), this outcome highlights the role of 


artihcial conduction in providing an effective description 
of the processes of mixing that lead to the redistribution 
of AGN feedback energy in central regions. 


In summary, we showed that our simulations produce a 
population of clusters that have a CC reasonable structure, 
both from a thermo- and chemo-dynamical point of view. 
Furthermore, once this result is achieved, the simulations 
also naturally produce a mix of CC and NCC clusters, 
which are similar to that observed. We trace the reason 
for the success of our simulations to the combined action 
of artihcial conduction, introduced to overcome numeri¬ 
cal limitations of standard SPH, and of AGN feedback, 
introduced to regulate star formation in massive halos. 

As a concluding word of caution, we would like to em¬ 
phasize that the success of the simulations presented here 
does not imply that we are self-consistently describing the 
details of all of the physical processes that lead to the cre- 




Fig. 3.— The effect of AGN feedback on entropy (upper panel) 
and iron abundance (lower panel) profiles for z = 0 clusters. In 
each panel the grey lines are the profiles of all clusters simulated 
without AGN. For comparison, we report the shaded area of the 
GC entropy (from Fig. 1, top-left panel) and the shaded area of the 
NCG metallicity (from Fig. 1, bottom- right panel). 



^/^500 

Fig. 4.— Effect of artificial conduction on the entropy profiles of 
the cluster shown in the upper panel of Fig. I^at two instances: z = 
0.8 (magenta curves) and z = 0 (blue curves). The simulation has 
been carried out with and without artificial conduction (solid and 
dashed lines, respectively), while leaving all the other parameters 
unchanged. The cluster simulated with artificial conduction was 
classified NCG and CC at the two epochs. 
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ation of cool cores: inflation of bubbles of high-entropy 
gas from the shocks of sub-relativistic jets launched by 
the AGN, gas circulation and turbulence triggered by the 
buoyancy of these bubbles, possible effects of magnetic 
fields in stabilizing them, thermal conduction, and cosmic 
rays. These are all effects for which we have circumstan¬ 
tial evidence from observations, but that are not included 
in our simulations. Still, our results demonstrate that our 
description of AGN feedback and artificial conduction pro¬ 
vides a realistic effective description of the above physical 
processes, without requiring any specific hne tuning. 
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